Single polymer dynamics: coil-stretch transition in a random flow. 
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By quantitative studies of statistics of polymer stretching in a random flow and of a flow field 
we demonstrate that the stretching of polymer molecules in a 3D random flow occurs rather sharply 
via the coil-stretch transition at the value of the criterion close to theoretically predicted. 



PACS numbers: 23.23.+X, 56.65.Dy 

The dynamics and conformations of an isolated flexi- 
ble polymer molecule in various flows form the basis for 
our understanding of hydrodynamics of dilute polymer 
solutions and remain an outstanding problem in poly- 
mer physics for several decades [1] . For the last 30 years 
these studies were concentrated on a well-characterized 
elongation flow by using classical techniques such as rhe- 
ology, light scattering, and birefringence [1]. In the re- 
cent years real-time imaging fluorescence microscopy of 
a single DNA molecule in a solution was developed to 
investigste its dynamics in elongation [2] and shear [3] 
flows, and also in their linear combination [4]. A good 
agreement with theory in all these cases was found. 

In contrast to it, studies of the polymer stretching in a 
general 3D random flow relevant to many physical and in- 
dustrial applications, and particularly to turbulent drag 
reduction, one of the basic problems of polymer physics 
[1], almost do not exist. The first experimental study of 
this problem was conducted only recently by macroscopic 
mechanical measurements [5] . It was shown that as a re- 
sult of a secondary random 3D flow superimposed on a 
primary applied shear flow between two disks, the poly- 
mer contribution to the shear stress increases as much 
as 170 times. If one assumes a linear elasticity of the 
flow-stretched polyacrylamide (PAAm) molecules, then 
the elastic stress causes 13 times polymer extension [5]. 
The main difficulty to perform an experiment to study 
polymer stretching and dynamics on a level of a single 
molecule is to create a random flow in a microscopic size 
volume. A conventional way to create chaotic or tur- 
bulent flow is to work at high velocities, V, and in a 
large size vessel, L, to reach large Reynolds numbers, 
Re — VL/v, where v is the kinematic viscosity. Large 
Re flow cannot be achieved in a hundred micron cell. A 
way to overcome this problem is to exploit a recently dis- 
covered elastic turbulence [6] . It is an apparently random 
flow that can arise in polymer solutions. 

Solutions of flexible polymers are visco-elastic fluids 
[1]. Polymer molecules stretched in a flow produce elas- 
tic stresses that grow nonlinearly with a flow rate and 
are manifested of many dramatic effects such as several 



orders in magnitude increase in extensional viscosity at 
high extension rates, the rod climbing effect, and partic- 
ularly the turbulent drag reduction [1]. The strength of 
non-linear effects is defined by the Weissenberg number 
Wi = tV/L. For elasticity-driven flows Wi plays a role 
analogous to that of Re for inertial turbulence in Newto- 
nian fluids. At sufficiently large Wi and arbitrary small 
Re a transition from laminar flow to the elastic turbu- 
lence occurs in flows with curvilinear trajectories [6]. The 
elastic turbulence bears the features similar to hydrody- 
namic turbulence [6]. The fluid motion is excited in a 
wide range of spatial and temporal scales, whereas sig- 
nificant increase in a flow resistance and effective mixing 
are observed [6,7]. Since the elastic stresses are indepen- 
dent of Re, they remain unaltered, when the vessel size 
decreases at the same Wi. Indeed, the elastic turbulence 
was recently observed in a curvilinear micro-channel flow 
[8]. Thus the elasticity-driven turbulent flow is a very 
appropriate system to study single polymer dynamics in 
a 3D random flow in a micro-device. 

In this Letter we report experimental observations of 
dynamics and conformations of a single DNA molecule in 
a 3D random flow and of the coil-stretch transition that 
accompanied it. 

The experiments were carried out in a swirling flow 
created in a d « 300/im gap between the edge of a uni- 
formly rotating glass rod of radius r\ = 1 mm and 0.14 
mm cover slip. The cell sidewalls were made of delrin 
with inner radius of r2 = 6 mm. The glass rod was glued 
into a metal shaft, which was driven via a belt by an 
optically encoded dc mini-motor with less than 1% rms 
velocity variation. 

We used two polymer solutions: One solution 
contained 80 ppm of high molecular weight PAAm 
(18,000,000; Polysciences) , and another had 36 ppm of 
A-DNA (48,502 bp)(Gibco). Both solutions have viscos- 
ity of about j] = 0.14 Pa-s at the shear rate 7=1 s _1 and 
the working temperature of 22.5 °C (Fig.lA). Two differ- 
ent solutions were used to generate the elastic turbulence 
and to demonstrate that the dynamics and conformations 
of a stained DNA molecule do not depend on the way a 
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random flow was produced. 

To study the dynamics and statistics of polymer 
molecules in a flow 10~ 3 ppm of A-DNA molecules, fluo- 
rescently labelled with YOYO-1 (Molecular Probes) at a 
dye/base ratio of 1:5 for > 1 hour, were added into both 
solutions. . Experiments were performed in a pH8 buffer 
consisting of lOmM tris-HCl, 2 mM EDTA, 10 mM NaCl, 
4%/3-mercaptoethanol, glucose oxidase (~ 50 fig /ml) and 
catalase (~ 10 fig /ml), 62% (w/w) sucrose (see e.g. Ref. 
[2]). At equilibrium the coiled A-DNA has R g — 0.73^m, 
while the entire contour length is I « 21/zm. So it may 
be considered as a "flexible" polymer with roughly 300 
persistence lengths [9,10]. Fluorcsccntly labelled DNA 
molecules were monitored from 15 to 100 /im above 
the cover slip via x63, 1.4NA oil immersion objective 
(Zeiss) with 0.4 //m depth of focus mounted in a home- 
made inverted epi-fluorescent microscope. Images of the 
molecules were digitized, and their maximum extension, 
R, was automatically measured. Various molecule con- 
formations in the turbulent flow that changed continu- 
ously were observed (Fig. IB). 

The data on the molecule extension were recorded at 
angular velocities ranging from fl = up to 7.0 s _1 for 
both cases of a laminar shear flow in a pure solvent and 
of a laminar shear flow via the elastic transition and in 
the elastic turbulence regime for polymer solutions. The 
results of velocity measurements as a function of time in 
the laminar flow as well as in the turbulent regime were 
shown in Fig.lC. 

Since we are interested in statistics of molecule stretch- 
ing, we present the probability distribution function 
(PDF) of molecule extensions, P(Ri), obtained by bin- 
ning the data for about 1000 molecules with 1 fim bins 
and normalized by the total number of molecules for each 
value of Wi in the both shear and turbulent flows (Fig. 2). 
The data on PDFs for 9 values of local Wi — T^lti/d 
at the local value of for the turbulent flow and for 
6 values of Wi for the laminar flow were obtained af- 
ter waiting initial transient period of more than 50 units 
of strain (7 = jt > 50) since change in Q was made. 
The smallest measured extension is ~ 1/zm that is close 
to R g . As Wi increases the shape of PDFs in the 
both laminar and turbulent flows changes dramatically 
but differently. In the shear flow PDF changes from 
strongly skewed at low Wi to rather flat and symmet- 
ric at high Wi. In contrast to this, in the random flow 
PDF varies with Wi from strongly skewed toward lower 
extension values to the skewed toward higher values. 
The scaled root mean square(rms) value of extension, 

° = [Ei(Ri-») 2 P(Ri)} 1/2 /» (where fi = £\ R 1 P(R l ) 
is the mean extension), remains constant on the level 
of ~ 40% in the pure shear flow, and reduced down to 
~ 30% in the turbulent flow as Wi increases. The scaled 
PDF skewness, s = [£;(#; - fi) 3 P(R t )}/a 3 , differs dra- 
matically for two flows: for the shear flow it reaches 



about zero at high Wi, while for the turbulent flow it 
changes from large positive to large negative values as 
Wi increases (see inset in Fig. 3). The latter reflects the 
fact that the statistics of polymer stretching change crit- 
ically, so that DNA molecules change from preferably 
non-stretched to mostly stretched as Wi increases. 

We also calculated the mean fractional molecule exten- 
sion, fi/l, as a function of Wi for both flows (Fig. 3). In 
the shear flow(data 5) fi/l increases gradually and reaches 
value of 0.44 at the highest Wi. Our data agree rather 
well with the data of Ref. [3] (data 2). In the turbulent 
flow fi/l rises sharper to the value of ~ 0.68(data 3,4)that 
lies between those reached in the shear(data 2,5) and 
elongation(data 1) flows (Fig. 3). On the other hand, the 
molecule extension in the elastic turbulent regime corre- 
sponding to the maximum of PDFs(data 6) shows much 
higher values, up to 0.85 (Fig. 3). 

Let us compare the results on additional polymer 
stretching in a turbulent flow obtained by the direct mea- 
surements and by estimations based on our early mechan- 
ical measurements for polymer molecules of different stiff- 
ness (DNA versus PAAm) [5]. From the single A-DNA 
molecule measurements it follows that an increase in \xjl 
due to a 3D random flow compared with the shear flow 
is about 4.5 times (Fig. 3). The estimates based on a lin- 
ear elasticity approximation and the measurements of the 
PAAm contribution to a shear stress in both shear and 
random flows give the factor 13 [5]. The discrepancy indi- 
cates, first, that the different stiffness results in different 
degree of stretching due to a turbulent flow, and second, 
that the assumption of the linear elasticity at such large 
extensions is not valid, and a nonlinear correction for a 
finite molecule extension should be taken into account. 
On the other hand, by using the worm-like chain model 
[11] that describes rather accurately force-extension re- 
lation for a DNA molecule, and using the extension dis- 
tributions in both flows one can estimate the change in 
the polymer contribution to the shear stress due to ad- 
ditional stretching in the 3D random flow. It gives the 
factor of about 60 that should be compared with the fac- 
tor 170 measured directly in the table- top experiment 
with PAAm molecules [5]. 

The main goal of the experiment was to verify whether 
a sharp variation in a molecular extension in a turbulent 
flow indicates the coil-stretch transition. More than 30 
years ago Lumley [12] first suggested a qualitative theory 
of polymer stretching in a random flow. It was recently 
revised, and the quantitative theory of the coil-stretch 
transition of a polymer molecule in a 3D random flow 
was developed [13,14]. Dynamics of a polymer molecule 
are sensitive to a fluid motion at the dissipation scale, 
where the velocity field is spatially smooth and random 
in time [12]. On this scale polymer stretching is deter- 
mined only by the velocity gradient tensor, k, that varies 
randomly in time and space. In a 3D random flow k 
always has an eigenvalue with a positive real part, so 
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that there exists a direction with a pure elongation flow 
[15]. The direction and the rate of the elongation flow 
change randomly, as a fluid clement rotates and moves 
along the Lagrangian trajectory. If k remains correlated 
within finite time intervals, the overall statistically aver- 
aged stretching of the fluid element will increase exponen- 
tially fast in time. The rate of the stretching is defined 
by the maximal Lyapunov exponent, A, of a turbulent 
flow, which is the average logarithmic rate of separation 
of two initially close trajectories. 

Stretching of a polymer molecule follows a deformation 
of a surrounding fluid element. So the statistics of poly- 
mer stretching in a random smooth flow depends criti- 
cally on A, and therefore, on the value of A • t. According 
to the theory [12,13] the polymer molecules should be- 
come vastly stretched, if the condition A-r > 1 is fulfilled, 
and the coil-stretch transition is defined by the relation 
X cr ■ t = 1 similar to that in a stationary elongation flow 
with the extension rate equal to A [16]. 

According to the recent theory [13] the tail of PDF 
of molecular extensions is described by the power law 
P(Ri) ~ R^ a , where a ~ (r _1 — A) in the vicinity 
of the transition. Positive a corresponds to the major- 
ity of the polymer molecules being non-stretched. On 
the contrary, at a < the majority of the molecules is 
strongly stretched, and their finite size is defined by the 
feedback reaction of the polymers on the flow [13] and 
by non-linearity of molecular elasticity [14]. Thus, the 
condition a = can be interpreted as the criterion for 
the coil-stretch transition in turbulent flows [13]. 

Solid lines in Fig. 2 are algebraic fits to the PDF tails. 
The exponents a obtained from the fits versus f2 _1 are 
presented in Fig. 4. The linear fit to the data intersects 
with a = at Q, cr — 0.68s -1 . 

In a separate experiment we measured statistics of A as 
a function of Q, around the elastic instability transition 
from laminar to turbulent flows. The measurements were 
conducted in the same set-up (besides the objective x20 
with 1.94 jum depth of focus was used) and in the same 
solutions but instead of the labelled DNA molecules ei- 
ther 1 /Ltm or 5.7 fim green fluorescent beads were added. 
A separation of two adjacent beads along Lagrangian tra- 
jectories was measured on many bead pairs (about 3000 
pairs at each value of fl, and each point represents the 
data collection during 30 minutes) as a function of O. 
From PDF of finite-time Lyapunov exponents, defined 
from the beads separation rate, the Cramer rate functions 
5(A) cx t^ 1 lnP(X,t i ) for several observation times U are 
found to collapse on one curve (see inset in Fig. 5) [17]. Its 
minimum defines the average Lyapunov exponent, A, for 
a given Q, [17]. The dependence of A on fl is presented in 
Fig. 5. From this plot one can clearly see that the elastic 
instability takes place at Q e i « 0.46s" 1 , while for the coil- 
stretch transition at Q, cr one finds X cr = 0.07 ± 0.017s _1 . 
It provides for the criterion of the coil-stretch transition 
X cr • t = 0.77 ± 0.20 that is rather close to the theoret- 



ically predicted value and the recent numerical simula- 
tions [17,18]. Small discrepancy can be attributed to the 
non-linear elasticity of a DNA molecule [14] and inhomo- 
geneity of the velocity field [13]. 
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FIG. 1. A. Schematic diagram of the set-up. B. Images of 
individual polymer molecules in a turbulent flow at 57 = 7s -1 
and Wi = 56.5. The horizontal bar indicates 5/i m. C. Tem- 
poral behavior of azimuthal flow velocity at two values of O: 
(1) 0.5s- 1 ; (2) 5.5s- 1 . 

FIG. 2. PDFs of polymer extension in turbulent (a,b,c,d) 
and steady shear (e,f,g,h) flows at several values of Wi. 



FIG. 4. a as a function of fl -1 . Black dots present the 
data after the coil-stretch transition, and open circles -before 
it. Solid line is the fit to the data, and dotted lines show the 
transition at a — 0. 



FIG. 5. Average Lyapunov exponent A as a function of 
n. The data are taken in: (I) PAAm and (2) A-DNA so- 
lutions. Inset: The Cramer rate function S(X) versus A at 
n = 0.83s -1 . 



FIG. 3. The mean fractional molecule extension calculated 
from PDFs as a function of Wi for various flows (explanation 
in text). Lines are to guide the eye. Inset: rms and skewness 
of PDFs as a function of Wi for turbulent (dots) and shear 
(squares) flows. 
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